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Abstract 

We present a set of three PDEs for three real scalars that are equivalent to the 
full Einstein equations without any symmetry assumptions. The main variables 
in this formulation are null surfaces and a conformal factor. 

Furthermore, for asymptotically flat spacetimes the free data (representing 
gravitational radiation) enters as the source term in the resulting equations. This 
could be important for an asymptotic quantization procedure. 

1 Introduction 

The Null-Surface Formulation of general relativity was developed in the 80's[T] and 
90'sj2l[3] with a double purpose. On one hand it provided a generalization of self-dual 
General Relativity (GR) |4J to full GR. It also presented GR not as a field theory but 
as a theory of surfaces on a four-manifold coupled to a scalar function. These surfaces 
were characteristic surfaces of some conformal metric. The scalar function played the 
role of a conformal factor converting the conformal metric into a vacuum metric (with 
a straightforward generalization if matter were present). 

The resulting field equations for NSF were given on the bundle of null directions, 
with the spacetime as the base space and the fiber given by the sphere of null directions. 
Furthermore, introducing the concept of spin weight in the fiber space the resulting 
equations were labeled according to this number. We recall that a sw-0 object is 
a scalar under rotations on the sphere, a sw-1 object represents a vector, etc. The 
original equations were then given as one sw-0 equation, two sw-1 equations and and 
two sw-3 equations. This set of equations were equivalent to the Einstein equations 
albeit in a completely different form and for completely different variables. 



1 



The purpose of this note is to show that one can reformulate NSF as a set of three 
sw-0 equations. It is remarkable that full GR without any symmetry assumptions can 
be written as three PDEs for three real functions on the bundle of null directions. 

In Section II we first give a brief review of the NSF including some results that 
differ the original derivations. However, most of the details are left for the references 
contained therein. We then introduce the new set of equations equivalent to the old 
ones. The details of the calculations are given in the appendix. In Section III we present 
peeling in NSF, i.e., the behavior of our variables as one approaches null infinity and 
show how the free data representing the gravitational radiation enters as the source 
term in the equations. We also show that the field equations adopt an extremely simple 
form if we assume that the spacetime points are close to null infinity and we give a 
perturbation procedure to obtain the null surfaces good up to second order. Finally 
we close this work with some comments on possible use of these results. 

2 A revised NSF. 

The null-surface approach to general relativity is a formulation based on the introduc- 
tion of two real functions, Z and Q, living on the bundle of null directions over M; i.e., 
onMxS 12 with x a in M, and (C, () on S 2 . These two functions capture the information 
contained in an Einstein metric: Z(x a , £, C) encodes the conformal structure of the Ein- 
stein space-time and singles out a preferred member of the conformal class of metrics, 
while Q(x a , (, C) represents the appropriate conformal factor that turns the preferred 
member into an Einstein metric. In the following, we show how these two functions 
are introduced. First, we present the original version and then the spin weight zero 
version of the equations, showing how the number of expressions is diminished from 
four complex equations to two real ones. 

2.1 The Kinematical Structure 

On a manifold M we consider an S^-family of functions Z(x a , C, C) an d its level surfaces 
Z = const, such that the equation 

g ab (x d )d a Z(x d , C, ()d b Z(x d , C, C) = (1) 

can be solved for a metric g ab (x a ) for all values of (CC)- The surfaces Z = const, are 
then null surfaces of the metric g ab , i.e., for every fixed value of (CC); the equation 
Z(x a X, — u represents a null foliation of the spacetime (M,g ab ). It is clear that 
kinematical conditions must be imposed on Z if we want a solution g ab of Eq. ([I]). To 
obtain the solution as well as the conditions we consider eq. ([T]) and small deviations 
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(C + d(, C + d(),of the equation, i.e., we take derivatives of ([I]) with respect to (C)Cj- 
By repeated differentiation with respect to £ and £, and using the independence of 
the metric on the variables (CiC) ; we obtain our two main results: the explicit metric 
components and the metricity conditions imposed on Z for a solution to Eq.([T]) to exist. 

The metric components are obtained on a natural class of coordinate systems (one 
for each value of (C? C) ) • This special class of null coordinates 9 l (x a , (, Q, i = 0, +, — ,1, 
is defined by the equation 

6\x a XX) = (Z,(5Z,BZ,mZ). (2) 

In the above equation, the operators S and S are the C covariant derivatives on 
the sphere holding constant x a [Sj. Eq. (P2J) should be interpreted as a coordinate 
transformation x a — > 9 l for every fixed value of (CC)- in the following we will use 
indices a for the ((,() independent coordinates x a and i for the (x a ,(,() dependent 
9 l . The inverse transformation is then given by 

x a = x a (0\C,C). 

Using the gradient basis d a 9 % ', together with the dual basis e" = diX a the metric 
g ab can be expressed in the new coordinates as g ab (x) = g^efej where the metric 
components g l i are explicitly given below. 

The g l i are all expressible in terms of two quantities, S 2 Z and g . By assumption 
g 01 is a positive definite scalar and can be written as g 01 = Q 2 . However, in the 
derivations it is more useful to keep its original definition, 

g m = g ab ( x a )d a Zd b mZ. (3) 

The function d 2 Z defines the complex scalar A via 

A(0\C,C) = &z(x a (0\,CC).X:O- 

As we will see below the metric components depend on <9jA and 3 and 9 derivatives 

The same operator S can be written in the 6 l coordinates in the following manner 

g = 5' + m'di (4) 
where 9' is the usual ( derivative holding 6 % fixed and is given by 

m° = e + 

m + = a (5) 
m- = e 1 

Q0 1 = g A _ 29+ . 
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Explicitly we have g lJ = g h lJ [A] with h tJ given by 



h 00 = h 0+ = h°- = o 

h m = -h + - = 1 
h ++ = h— = -d x A 
l I A- §IdiA 



where 



h +1 = h 



dxS A + (h ij - h^SlSDdikdjk + br*d&k + h +i d { dA + 2 



2 (1 - i^A^A) 



(6) 



(1 + i^A^A) 



A = -d^A + a + A - h—d_A. 
The metric can then be written as 



(7) 



2.2 The metricity conditions. 

The scalars g 01 and Z are not arbitrary nor independent of each other. They must 
satisfy two coupled differential equations, referred to as the metricity conditions. The 
first one is obtained applying the 9 operator to g 01 , i. e., to the expression (j3J 

&y 1 = g ab {d a e + d b e 1 + d a e°d b M) (8) 

which in components gives 

cV 1 = g 01 (dJiA + h +1 ) . (9) 
Inserting explicitly the value of h +l in Eq.(j9]) gives 

dj5A + d + A = Mi (10) 

with 

M I = 2\1- -diAdiA J 9 In # 01 - S^Ad-A - -(SiSA - 3_A - d x Ad+A - d{Rd{dh)dih. 
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For the second one we apply the <3 3 operator to the expression (JT]) obtaining, 

g ab (3d a 9 + d b A + d a 9°d b M) = (11) 

which in components gives 

d^A - 3d_A = M n . (12) 

with 

M n = 3diA (<9+A - h +1 ) . 

The two metricity conditions constitute the requirement on Z in order for a metric 
g ab to exist and not depend on (£, £) and such that Z = const are characteristic surfaces 
of the metric, for every value of (C>C)- They leave g 01 undetermined up to a factor 
dependent only on x a . The second condition has another geometrical meaning and it 
is related to seminal work done by E. Cartan on solutions of certain PDEs. To briefly 
review this second point of view we start our construction from the equation 

3 2 z = A(z,dz,Bz,mz,(,0 (13) 

where A is assumed to be a given function of the six coordinates (9 l , £, Q which satisfies 
the integrability condition 

9 2 A = 9 2 A, (14) 

Note that the coordinates x a have disappeared. If we analyze the solution space of this 
equation we find that there are four constants of integration and the solution can be 
written as 

u = Z(x a ,(,() (15) 

(the spacetime points reenter as the solution space of this master equation). Further- 
more, if we want to characterize the special class of functions that are equivalent to 
each other via diffeomorphisms u' = F («, C>C) an d C' = ^ (C> C) (called fiber pre- 
serving transformations) we find that they should give the same conformal structure 
on the solution space, i.e., they should give the same conformal metric defined above. 
The last condition f lT2|) is then the requirement that the normal metric connection has 
vanishing torsion [B]. 

From our point of view, this fiber preserving transformation Z' = F (Z, £,G and 
C = G (£, C) clearly gives the same conformal structure since for a fixed value of (£, C) , 
Z = const, and Z' = const, are both null surfaces for the same metric. 
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2.3 The real (sw-O) metricity conditions. 

We note that equation fl9]) is a complex pde for a single real variable g 01 . In principle 
there could be no solutions unless the integrability conditions are identically satisfied. 
This is indeed the case. Taking 3 of Eq. (jSJ), we obtain: 

28CV 1 = g ab (d a Q 2 Ad b 6° -d a Ad b A), (16) 

or 

m g 01 =^g°\d 1 d 2 A-ti j d i Ad j A). (17) 

Thus, B<Sg 01 = 59g 01 if 9 2 A = 5 2 A, but this is true from the starting assumption, 
equation (JHJ). 

Note also that we could use the real Eq. (fT7l) instead of (jSJ) if we are interested in 
regular solutions for g 01 . This follows from the fact that the only regular solution of 
5F 1 = for a s.w.l function F 1 is F 1 = 0. Thus, we could either use or (I17p as our 
first metricity condition. 

With this in mind, we expect the same feature from our second metricity condition, 
Eq. ( I12p . Directly from 

m 3 g ab (x d )d a Z(x d , C, Od b Z{x d , C, C) = 

(see the Appendix I for detailed calculations) we get 

9i858 2 A + 6d x 8 2 A + 9h ll dS 2 A + 3/i + ^5 3 A + 3/r fc <9 fc o :3 A+ (18) 
+ h ij (QdiAdjA + 9dMdjBA + d&KdfiA + Zd&BAdjK + W^AdjA) = 

Again, and using the same argument as before, if we are interested in regular 
solutions for A we can take Eq. ffTHj) instead Eq. f lT2"j) . 

Summarizing, given two functions A, g 01 that are solutions to ([IT]) and (FlBl . the 
level surfaces Z (x a , C> C) = const, of the solution to (TT3T) are null hypersurfaces of the 
metric defined by Eq. ([1]). In the following sections of this work we will use the spin 
weight zero real formulation, equations (1171) and (1181) to derive the field equations for 
Z. 

2.4 The field equations 

On this kinematical scheme we impose the (trace-free) vacuum Einstein equations by 
Rij — \gijR = 0. These are a set of nine equations for the metric components g^j in 
our coordinate system 9 % {x a , £, C) for arbitrary but fixed (C, C)- However, we can show 
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that these equations are generated by just one equation in the null bundle coordinates 
(#\ C> 0- Consider the following contraction 

R ab (x a )d a Z(x d , C, ()d b Z(x d , C, C) = 0. (19) 

Taking 5 and 6 on this equation we generate the remaining eight trace-free vacuum 
Einstein equations in an analogous fashion as was done on the previous subsection. 
Writing f|T9|) in terms of A and g 01 = Q 2 gives: 

dfn = ORu(A) 

with i?n(A), a component of the Ricci tensor of h % \ given by 

R n (A) = R ab (A)d a Z(x d , C, Qd b Z(x d , (, () = l^A^A + ^(d lQ ) 2 - ]-d 2 q, (20) 



q = 1 - d 1 Ad 1 A (21) 
The full set of Einstein equations can then be grouped together as 

<9^-ORh(A) = (22) 

2§g(fi 2 ) - tf^d^A - WdiAdjA) = (23) 

^598^ + 6a!5 2 A + 9h u diB 2 A + 3h +j djB 3 A + 3h- k d k d 3 A+ (24) 



+ h ij (GdiAdjA + 9dMdjdA + dMdfiA + 3^g9Aa,A + 3dMAdjA) = 

The full system consists of two coupled PDEs for Q and A plus the equation ( )24|) for 
A. From now on this set of equations will be referred as the NSF equations. 

Note that the first NSF equation, one real equation for (Q,A), is equivalent to the 
trace-free, Ricci flat equations, nine second order equations for six real variables in 
the standard formulation (plus identities). However, the metricity conditions must be 
incorporated in the NSF since they give the conditions that the level surfaces Z = const. 
are indeed null surfaces for the metric of the spacetime. Thus, the three real equations 
of the NSF are equivalent to the equations that yield Einstein spaces in GR. 

3 Asymptotic NSF 

Though any vacuum Einstein space-time can be investigated in this manner we make 
the specialization, here, to asymptotically flat vacuum space-times. In this case the 
geometrical meanings to the various quantities become more focused and clearer and 
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the differential equations become easier to handle. We begin by assuming that the 
spacetime is asymptotically flat along null directions with a null boundary X + . Null 
infinity can be coordinatized with a Bondi coordinate system, (u, C, C) with u £ R, the 
Bondi retarded time, and (£, C) £ S 2 labeling the null generators of X + . With this 
notation we can give a precise meaning to the null surfaces described by Z(x a , £, £) = 
u = const; they are taken to be the past null cones of the points (u, (, Q of X + . If in 
addition we set (jZ(x a , £, Q = uj = const, we isolate a null geodesic in this light cone. 
(Any null geodesic can be labelled by these five parameters, the apex (u, (, Q of the cone 
at I + , and its "null angle", (u,u).) Furthermore, setting dBZ(x a X,Q — R — const. 
we describe a point on this null geodesic. 

In addition to this meaning of Z, there is a dual meaning, namely, if the space-time 
point x a is held constant but the ((, Q is varied over the S 2 , we obtain a two-surface 
on X + , the so-called light-cone cut of X + . It consists of all points of X + reached by 
null-geodesics from x a . Z is then referred to as the light-cone cut function. We now 
another interpretation, not only of Z{x a , (, £), but also of both, oj = <5Z(x a , £, () and 
R = (5dZ(x a ,(,C)- w is the "stereographic angle" that the light-cone cuts make with 
the Bondi u = const cuts and R is a measure of the curvature of the cut. 

3.1 Peeling in NSF. 

Consider a Bondi slicing of X + with the associated coordinates (u, £, C) and asymptotic 
shear o~b(u, (,() for the associated Bondi null surfaces; then the asymptotic shear o~z 
of a light cone with apex at x a evaluated at X + is given by 

o- z (x a XX)=o- B (Z,(X)-d 2 Z , (25) 

where Z (x a ,C,C) is the light cone cut from x a . This relationship is known as Sachs' 
theorem. From equations (ITHjl and ( 125]) we see that 

A = a B (u, C,C) - cr z (u,u,u, R,(,0, 

where we have used the coordinates (u, oj, uj, R) to label the apex x a . We want to study 
the behavior of A as we move the apex x a along a fixed null geodesic characterized by 
a given value of (u, oj, uj, £, Q, i.e., along R. Since our spacetime is asymptotically flat 
it is clear that as the apex R — > oo its light cone resembles more and more like a flat 
one and thus it should have a vanishing shear in this limit. Thus 

lim A = a B (u,C,Q. 

R^foo 

Actually, this intuitive argument can be made precise if one uses the optical equa- 
tions for the shear and divergence of the light cone from x a . One can show that (see 
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Appendix 2) as R — > oo the leading order term in oz is given by 

Z 12R 2 

and thus 



+ O (R- 3 ) , (26) 



A = ^ + I ^ + 0(i?- 3 ). (27) 

Remark: the fact that lim^oo A = cr B (u, (, Q shows that A explicitly depends on the 
free Bondi data. Moreover, Eq.( l2T)) can be used to obtain all the metric components 
and derived tensors in the asymptotic limit. 

Other scalars have interesting asymptotic behavior as for example SA or S A. To 
see this we write the 8 operator in the 6 % coordinates as 

B = tf + ua u + Rd u +~fcd CJ + (tiK-2u)d R . (28) 

Taking <3 of Eq. (J27D we get 

5A = da B + (R- 1 ) (29) 

2 

and 5 ( 1271) gives: 

B 2 A = ta B + O (R°) . (30) 

We can make no asserts about this last decay just from kinematical considerations. 
After we derive the field equations below we obtain the explicit form of the second 
term on the r.h.s. of the equation. 

3.2 Field equations and free data. 

To derive the field equations we start with Eq.( l23|) . i.e., 

d R d 2 A = 2^- +h ab d a Ad b A. (31) 

Integrating (|3T|) over R we obtain 

r°° Bar) 2 

B 2 A + j R [2— +h ab d a Ad b A]dR'=(d 2 A) R ^ oo , (32) 

where (B 2 A) R _ S>00 is a function that depends on (w, u, aJ, (, (). Using the remaining field 
equations plus the peeling behavior of all the scalars involved one can show that (see 
Appendix 3 for details) 
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(aJPA)^ = tfa B + <5 2 W B + o B o B . (33) 
Using the fact that d u commutes with 5 2 and 9 in the limit R — > oo we finally obtain 

/u 
a B W B du. (34) 
-oo 

The constant of integration has been set equal to zero since we restrict the BMS group 
to the Poincare group by asking a -)• as w -> -oo. The integral term cr B a B du 
is the total energy radiated by the system and it is finite by assumption^ 
The final equation reads 

8 2 A + / [2—^ +h ij d i Ad j A]dR' = 3 a B + Q 2 a B + / o B o B du (35) 

JR " J-oo 

where Q must be solved from the Einstein's equation f )22|) . Equation fl35|) together with 
( 124"|) and fT22|) is the full set of NSF equations for asymptotically flat spacetimes. Note 
that Bondi shear cr B (u, £, Q enters as a source term in Eq. fl3"5l) . This shows the explicit 
dependance of A on the free data. 



3.3 Leading Order Equations 

We have shown above that lim^oo A = <r B (u, (, Q. We now study the behavior of 5 2 A 
in this limit. For this we study the integral term on the l.h.s. of Eq. (l35l) . From Eq. 
( f20|) one can see that 

Q(A) = 0(R- S ) + ... 
It then follows from Eq. (J225 that 

Q = 1 + O (R- 6 ) (36) 

One can also see that 

WdiAdjA = O (R- 3 ) (37) 
Thus, the leading order term for 5 2 A on the region R » 1, is given by 

§ 2 A = 5 2 cr B + d 2 a B + / o B o B du + O (iT 2 ) . (38) 



1 Note that there exists a difference in the pius sign of the las term of Eq. (|34|) and the minus sign 
in the same term of the peeling of the main equation (15) in [7] 
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Note that if we keep the leading order term in the above equation and we replace u 
and A for Z and 3 2 Z respectively we get a field equation directly for Z, i.e., 

-2 r z 

S 2 c1 2 Z = 9 o B + &a B + / a B a B du. (39) 



Eq. (1391 is a real generalization of the good cut equation of if-space [I]. It also 
resembles the nice section equation of Moreschi [8] and the Geroch-Winicour equation 
P|. Note however that the above equation is different from the ones mentioned before. 
The final conclusion is that the solution space of Eq. f l3"§|) is different from the ones 
obtained in the other formulations since there is a one to one correspondence between 
points in the solution space and cuts of null infinity. An analogous statement would 
be that the light cone cuts produced by points of an asymptotically flat spacetime will 
satisfy Eq. (I31?l) in the limit R —¥ oo. 

3.4 A perturbative solution for A. 

Since the source term in Eq. fl33|) is quadratic in the Bondi data a B we want to obtain 
here an equation for which is correct up to this second order. To do this we star with 
Eq. (1221) and solve perturbately for Q, obtaining 



n = l + d R 'Ad R 'A + J J d R „,Ad RI ,,AdR!"dR" + 0(A 3 

R' R" 

Thus, the field equation for A reads 



S 2 A+ J [35 I d R ,Ad R ,A + J J d R ,„Ad R „,AdR'"dR" +r ] ij d l Ad j A}dR' = d 2 a+d 2 a+ j a B a B du 



R \ R' R" 



(40) 

where rf- 7 represents the Minkowski spacetime metric. This equation can be solved 
perturbately as follows, we first start with teh zero order solution for Z^°\ 

a® = 
A<°> = 
9 2 Z(°) = 

= x a e a (4i) 



n 



The first order term A* 1 ) is then obtained from 

8 2 A« = 8V(z(°\c,C) + S 2 a(^ (0) ,C,C) 

6 2 [A« - a(Z^ , C, C)] = d 2 a(Z^ 
and integrating at x a constant we get 

A« = a B (Z(°),C,C) + J Ki & 2 a(x a l' a ,CX)dS 2 (42) 

5-' 

where K 20 is the corresponding Green's function for the B 2 operator [5]. Finally 

a(Z^,U)a(Z^,COd 

-oo 

pR / POD POO 

+ / [m(d R A^d R A {1) ) +m( / d 2 RIII A^d 2 RIII A W dR'"dR 

Joe ^ ' ii?" 



(43) 



+ rf 3 d l A {l) d j A {l) }dR' 

As we can see, the expression (1431) is a PDE on the sphere for the cut function Z^> 
is completely know from previous integrations. The equation can be integrated using 
the corresponding Green's function for 8 2 5 2 . 



4 Summary and Conclusions 

In this work we have obtained a new derivation of the field equations for NSF. By 
proper differentiation we obtain spin weight zero equations that are equivalent to the 
original ones. The resulting equations constitute a set of three real PDEs for three 
real functions on a six dimensional space. The solution of those equations yield null 
surfaces on the spacetime, (equivalent to the knowledge of a conformal metric) and the 
conformal factor that yields a vacuum metric. 

For asymptotically flat spacetimes one obtains equations where the source term is 
the free Bondi data for gravitational radiation. This new result is particularly impor- 
tant for several applications. Since the light cone cut Z(x a ,C,C) °f an Y point x a of 
the spacetime must satisfy this equation, one can define a family of light cone cuts 
associated with a worldline by simply giving the parametric form of the worldline as 
x a (r) with r the proper time of the worldline. In particular, any candidate for a center 
of mass worldline should be a solution of Eq. ( 135]) . In this context, none of the ap- 
proaches that are present in the literature satisfy this equation and it should be useful 
to incorporate Eq. (135|) in those approaches. 
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It is also worth pursuing further studies on the solution of Eq. fl39|) . We have the 
machinery needed to investigate the dynamical properties of the spacetime defined on 
the solution space of this equation. 

Finally, we would like to discuss and issue that must be addressed each time a new 
variable is introduced to describe General Relativity. What is the advantage, if any, to 
present this non local variable instead of the metric of the spacetime? After all, from 
knowledge of the metric one can construct null geodesies, null cones, etc. 

So far one can point out some advantages that come directly from the NSF for- 
malism; Since the equations yield global variables any perturbation (or numerical) 
calculation of the solution modifies the null surfaces of the spacetime at each order of 
the perturbation calculation, a task that is virtually impossible with standard variables. 
Another advantage of the NSF is that the free data enters directly into the equations 
with a clear physical meaning, namely, the data is the gravitational radiation that 
reaches null infinity. One of the outstanding problems in a Cauchy formulation of GR 
is to define what constitutes physically relevant data and this issue is non existent in 
the NSF. At a classical level one can foresee the role of NSF not as competing with the 
standard approach but rather as an alternative formulation available to discuss global 
issues and/or solutions without any symmetry. 

The perception as to which is a good variable changes completely when dealing 
with a quantum theory of gravity. In Loop Quantum Gravity the metric field is not 
an observable since it is not a self adjoint operator. Instead, non local variables such 
as holonomies around loops play a leading role in the theory. Here NSF could play an 
important role for several reasons. By construction the main variables are non local 
observables and the free data Ob is promoted to a quantum operator erg that obeys free 
commutation relations. Since lim^oo A = C> C) one nas i n NSF a well defined 

variable that asymptotically approaches the free quantum operator a B . This could be 
particularly useful for an S matrix approach to quantum gravity. 

Appendix I 

In order to obtain equation ffTSl) we apply 6 3 <3 3 to g ab {x d )d a Zd b Z = which corre- 
sponds to applying 5 3 to ffTTj) 

d 3 g ab {3d a 9 + d b A + d a e°d b QA) = 0, 

obtaining, 

= g ab (d a d 3 dAd b u + 3d a d 3 Ad b u + 3d a e 2 clA<9 6 ZJ + 3d a mAd b A + 3d a mAd b A- 
- 18d a Md b uJ + 9d a WAd b R + <9 a 6A<9 6 c3A - Qd a Ad b A + 9<9 a gA<9 6 5A). (44) 
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Performing [5, 5] F s = 2sF s and inserting 5 2 A = 5 2 A, equation (|44p can be put in the 
form 

= g ab (d a Bm 2 Ad b u + 6d a d 2 Ad b u + 9d a B 2 Ad b R + Qd a Ad b A + 9d a BAd b M+ 
+ d a Md b dA + 3d a B 3 Ad b u + 3d a d 3 Ad b uj + 3<9 a 5SA<9 6 A + 3d a 83A<9 6 A) (45) 

which is clearly a real equation. Setting f l4^|) in components we have, 

= diBSo^A + 6di8 2 A + 9h u d$ 2 A + 3h +j djd 3 A + 3h- k d k & 3 A+ 
+ h ij (QdiAdjA + 9d$~AdjdA + dMdfiA + 3d i : &dAd ] r A + 3^55X5^) 

which is Eq. fll8p in the main text. 



Appendix II 

In this Appendix we calculate the form of az(u, cu, u, R, (, (). We start from the 
equations for the optical scalars associated with null geodesic congruences on a vacuum 
spacetime 

% = P 2 + °z°z ( 46 ) 

^ = 2pa z + ^ , (47) 

with S being the affine parameter and ipo = C1313 as defined in the Newman-Penrose 
formalism [TO] . 

In order to integrate the equations (H6l) and (T47j) we introduce a perturbative scheme 
in terms of a strength factor, considered as e< 1, such that we can write 

CT = (7 (o) +e(7 (i) + e V 2 ) + ... 

p = p(°) + e p« + e y 2 ) + ... 

Then, we solve through a perturbation procedure. 

Order zero: Solving equations (|46j) and (1471) for Minkowski space 0% = and 

^ = (p(0)) 2 

whose solution is 



p(°) = — 
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with So representing the position of the apex of the light cone. 
Order e: Writing equations (I46p and ( l4"T|) for and as 



~~dS~ 



(0)n(!) 



2p^p 



dS 



we find that the solutions are 



a 



(i) 



(s - s y 

Taking the limit lim^oo S 2 a^ , 



p (1) = 







1 


2S 


C'2 


2S 2 + 


3S 3 


AS 4 



+ 



The last step is to show that in the limit lim^oo ^ — > 1. This is true by virtue of the 
defining equation for g 01 , namely 



dR 

Is 



9 



01 



asymptotic simplicity implies lims^oo g 01 — > 1. Thus, the behavior of &z is 0(R~ 2 ), 
where R is the position of the apex. This last expression corresponds to equation ( 1261) 
of the main text. 



Appendix III 

In this Appendix we obtain equation ( J33l) . Although the intervening equations are 
long and the derivation cumbersome the only goal is to obtain a limit using the peeling 
behavior of our variables. Most of the time this peeling is easy to take into account 
but there is one term that, as we see below, requires a bit of effort. 

We first split the linear and upper order terms in A in Eq. (T2"4"j) as, 

di8S8 2 A + 6di8 2 A + 9d u d 2 A - 3ob8 3 A - 3d w 5 3 A+ 

+ 9(/i 1+ <9 + 5 2 A + h^d^A + h 11 d 1 d 2 A)+ 

+ 3{h ++ d + B 3 A + /i +1 di8 3 A + h— <9_3 3 A + h~ 1 d 1 d 3 A)+ 

+ h ab (6d a Ad b A + 9d a dAd b 8A + d a Md b dA + 3d a dMd b A + 3d a mAd b A) = 0. (48) 



15 



We use the commutation relations between ft and 8 written as 

[ft, 8] = 6+ {do - 2ft) + Sid- + dM + + dMdi = [ft, 5]° + [ft, 3] 1 

and its complex conjugate, where [ft, 5]° = <^ + (ft) — 2ft) + 8}d~ and [ft, 3] 1 = ftA5 + + 
ftSAft to reexpress the first term of (l4"8j) . ft838 2 A as 

ft838 2 A = 83ft8 2 A + d_8 3 A + d + 38 2 A - ft,8 2 A + 2ft8 2 A + 
+ [ft , 8] 1 38 2 A + 8 [ft , 3] 1 8 2 A + [8, dJ\ 1 8 2 A 

and replacing in (T4"8l) we get 

4d u 8 2 A - ft^ 3 A - ft,3 3 A + 4a R 8 2 A + ^Sft^A + ^M IIRe = (49) 

with 

M IIRe = 9 (h 1+ 8 + B 2 A + h^dJ&A + /i n ft3 2 A) + 

+ 3 (h ++ d + B 3 A + h +1 d x WA + /zr-a_5 3 A + /^ftS 3 ^) + 

+ h ab (6d a Ad b A + 93 a 3Aft8A + <9 a 8Aft3A + 3<ft38AftA + 3<9 a 39A<9 6 A) + C(A). 

(50) 

where 

C= [ft,8] 1 88 2 A + 8[ft,8] 1 8 2 A + \E,d„] l d 2 A (51) 
= [ft , 8] 1 3 3 A + [ft , 3] 1 8 3 A + [8, <9_] 1 8 2 A + [8, 8+} 1 3 2 A+ (52) 
+ ( [8, ft] 1 A8+ + [8, ft] 1 3Aft + ftA [8, d+] 1 + ft3A [8, ft] ^ 8 2 A (53) 

= [ft , 8] 1 3 3 A + [ft , 8] 1 8 3 A + [8, 9_] 1 8 2 A + [8, 3 2 A- (54) 
- [(ftA<9_A + ftSAftA) a + 8 2 A + (ftA<9_8~A + ftSAftSA) ft3 2 A + cc] . (55) 

From the equations ( )22l) . ( 1271) . ( 13TI) . and ( 1501) . we can easily determine the decay 
in R of the three last terms in ( 1491) . This is not the case for ftj8 3 A + ft,3 3 A and an 
auxiliary expression must be found. It is possible to relate these two terms with the 
free data by a manipulation in the two metricity conditions. First, we equate the spin 
weight of the metricity conditions performing 

3cUHn} : 

3ftjft^3A + Sft^oU = 3ftjM/ 

d„m ■ 

d R d u M - 3GU.A = d u M n 
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adding these two, we get 

d R (3c^BA + d„M) = 3c^M/ + d u M u 
and integrating over the variable R, 

r R 

3d UJ QA + d ul QA = 4:a B + / (3flbM/ + d u M n ) dR'. (56) 

J oo 

We see that Eq. ( 15 6 p is a s.w.2 expression and applying 8 2 we turn it into a s.w.O 
one. If we carry out the corresponding permutations and the substitution of 8 2 A = S 2 A, 
Eq. (j56p can be reexpressed in the following manner: 

3oL9 3 A + d u &A - 63 U 8 2 A + !2d R d 2 A = 4S 2 a B + N (57) 

where N = -8~ 2 f* (35b-Mj + d u M n ) dR'-3 [d u , 8] 1 8A+ [d u , 8] 1 88A+8 [d u , 8] 1 8A+ 

38 [3*8] ^ + 3 [6br, B] 1 S 2 A + 6 [a i? ,8] 1 8A + 4 [^Sj'A. 
Adding to (1571) its complex conjugate expression we obtain, 

6b8 3 A + d^A - 3d u d 2 A - Qd R d 2 A = B 2 a B + d 2 a B + N Re 

with Njte = \ (N + N) representing second and upper order terms in A. Inserting this 
last expression in ( 149)) the following result is obtained 

d u B 2 A + d R &A + ^W R d 2 A = fa B + d 2 a B + N Re - ^M IIRe 

After a carefully study of the peeling of every term in this last expression we write 

(S u 8 2 A)^ oo = ta B + d 2 a B + a B a B 
which correspond to equation (1331) in the main text. 
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